#include <iostream>
#include "dnamap.h"
using namespace std;
using namespace dnamap;

int main()
{
    dnamap::dnaseq myseq;
    string xin, a, b, c;
    cin>>xin;
    myseq.enter(xin.c_str());
    cout<<myseq.nucleotide(0, 990)<<endl;
    protrans mytrans("codon.txt");
    a = myseq.getframe(0);
    cout<< mytrans.translate(a) << endl;
    restriction ecor("EcoRI", "GAATTC");
    sitelist mysite = mytrans.synmut(a, ecor.siteseq);
    dnamap::dnaseq mutseq(a);
    cout<<mutseq.nucleotide(0, 990)<<endl;
    return 0;
}
